Hubert's 17th problem and the quantumness of states 



in 
o 
o 

(N 

Ph. 

< 

O 
(N 

cn 
> 

OS 
(N 
O 

oo ; 

o . 

O ' 

Oh! 

> ■ 

a ■ 
^ . 



J. K. Korbicz^, J.I. Cirac^, Jan Wehr'^ and M. Lewenstein^ 
^ Institut fur Theoretische Physik, Universitdt Hannover, D-30167 Hannover, Germany 
Max-Planck Institut fur Quantenoptik, Hans-Kopfermann Str. 1, D-85748, Garching, Germany and 
^ Department of Mathematics, University of Arizona, 
617 N. Santa Rita Ave., Tucson, AZ 85721-0089, USA 

A state of a quantum systems can be regarded as classical [quantum) with respect to measurements 
of a set of canonical observables iff there exists (does not exist) a well defined, positive phase space 
distribution, the so called Galuber-Sudarshan P-representation. We derive a family of classicality 
criteria that require that averages of positive functions calculated using P-representation must be 
positive. For polynomial functions, these criteria are related to 17-th Hilbert's problem, and have 
physical meaning of generalized squeezing conditions; alternatively, they may be interpreted as non- 
classicahty witnesses. We show that every generic non-classical state can be detected by a polynomial 
that is a sum of squares of other polynomials (sos) . We introduce a very natural hierarchy of states 
regarding their degree of quantumness, which we relate to the minimal degree of a sos polynomial 
that detects them is introduced. Polynomial non-classicality witnesses can be directly measured. 



In recent years there has been a lot of interest in clas- 
sifying states of quantum systems with respect to their 
quantum nature. In particular, the problem of character- 
izing entangled states has attracted a lot of attention 1| 
because of its vital importance for quantum information 
processing. One of the aspect of this problem concerns 
non-locality of quantum mechanics and violations of Bell- 
like inequalities 2] and existence of local hidden variable 
models. The problem of existence of a classical proba- 
bilistic description of quantum states of a single system 
has, however, longer history and can traced back to the 
seminal papers of Glauber and Sudarshan 

Let us consider a harmonic oscillator Hilbert's space, 
and fix the canonical creation and annihilation operators. 
In the Refs. 3] it was shown that any state g can 



a, a 



be uniquely put into a form diagonal in coherent states 
|a): 



J d'^oeP{a, a*)\a){a\ , 



(1) 



where a = x -\- iy, (Pa = dxdy. The integration with 
P in 1^ is be understood in the distributional sense 0. 
Hermiticity and normalization of g imply that P* = P 
and ^ d^aP(a,a') = 1, while positivity implies that 
/ d^aP(a, a*)|(a|V'}|^ > for every?/'- 

A state of a quantum systems is classical with respect 
to measurements of a given set of canonical observables 
iff the Glauber-Sudarshan P-representation is a well de- 
fined, positive phase space distribution 5]. Mathemati- 
cally speaking, in such a situation P defines a probabilis- 
tic measure /z on the phase-space through: 



^ fi{n) / Pd'^a > 



(2) 



Statistical properties of state possessing the positive P- 
representation are as those of the classical statistical en- 
semble, described by the measure fi; that explains why 



such states are called classical. However, the class of al- 
lowed P's is larger than that P|, and there exist non- 
classical states (such as squeezed, or Fock states) for 
which the latter integral does not always exist, or attains 
negative values. 

In this Letter we derive a family of classicality crite- 
ria that require that averages of positive functions cal- 
culated using P-representation must be positive. For 
polynomial functions, these criteria are related to 17- 
th Hilbert's problem, which in its simplest form states 
that not every positive semidefinite polynomial must be 
a sum of squares (sos) of other polynomials Q . Our cri- 
teria have physical meaning of generalized squeezing con- 
ditions, and may be interpreted as non-classicality wit- 
nesses (in analogy to entanglement witnesses, |3]). We 
show that every generic non-classical state can be de- 
tected by a sos polynomial. 

Let us begin by observing that the set of probabilistic 
measures forms a convex subset of the set of all P's. The 
extreme points of this set are point concentrated mea- 
sures {5^{a — f3);f3 e C}, and the decomposition into 
these points is unique Q . Hence the classical states form 
a generalized simplex A. This is the general feature of 
sets of probabilities for classical systems 0. Therefore, 
geometrically the problem of distinguishing between clas- 
sical and nonclassical states amounts to the operational 
description of the simplex of positive measures A in the 
space of all P-distributions. We note that one encounters 
closely related problem in the study of quantum entangle- 
ment in multipartite systems (0 and references therein), 
with the difference that the convex subset of classically 
correlated states is not a simplex. 

The solution of such statedproblem was recently pro- 
posed by Richter and Vogel [l^. They studied the char- 
acteristic function of P, i.e. its Fourier transform: 



Ho 



tv{g:WiO :}, 



(3) 
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where ^ = ^r+if » and W(^) — e^° is the Weyl opera- 
tor. In what follows the hat," , will always denote Fourier 
transform. The criterion detecting positive measures is 
then provided by Bochner's theorem 

P is a Fourier transform of a probabilistic measure iff 
P is of positive type, i. e. for each number n and all 
possible sets fi,...,Cn G the n x n matrix P^j := 
P{^i — is positive semidefinite (psd). 

The further test of Py being psd for fixed ^i, . . . , ^„ is 
carried out using determinant criterion: a n x n matrix is 
psd iff determinants Dfc, k = 1 ... 71 of all of the principal 
submatrices are non-negative. This finally leads to the 
hierarchy of conditions: a state g is nonclassical iff there 
exist k > 2 (for k = 1, Di = 1 due to normalization) and 
points £_i, . . . such that Dk < 0. 

Our solution of the classicality-quantumness problem 
follows also from the Bochner's theorem. Note that the 
condition for P to be a function of positive type can be 
equivalently rewritten as: P is a function of positive type 
iff for aU X e P(M2), J d^ad"^ Pxia)* Pia - /3)x(/3) > 0, 
where X'(E^) is a space of smooth test functions with 
compact support |11|. Using the convolution theorem 
the last integral is equal to / d'^aP{a,a*)\x{—ct)\'^. From 
Fourier transform theory, x can be analytically continued 
to a function from Z(C^), the space of entire functions, 
satisfying specific bonds 'l2l|, and every element of Z(C^) 
is of that form 11, 12] . Hence we obtain the criterion for 
classicality [l^: 

P defines a probabilistic measure iff 



V/eZ(C2), d^aP\fMf>0, 



(4) 



where /r denotes the restriction of f to . 

Our approach offers new insights into the problem, and 
connects it to the methods used in the study of separa- 
bility. From 10} we obtain that a state g is nonclassi- 
cal iff there exists a test function / £ Z(C^) such that 
/ d^aP I/rP < 0. Since /r is real-analytic this condition 
can be rewritten as: 



tr{^): |/R(a,a^)|2 :} < , 



(5) 



implying that the state is nonclassical iff there exists an 
observable : \ f's,{a,a^)\^ : detecting it. Geometrically, the 
condition J d^aPj/Rp = defines a hyperplane in the 
set of all P distributions and hence a state is nonclassical 
iff there is a hyperplane separating it from the simplex 
A. This is essentially the same approach as the one used 
in the theory of entanglement witnesses Hi There- 
fore, we propose to call the observable from expression 
(Is} nonclassicality witness. The above approach can be 
generalized if we allow the test functions / to depend on 
the state g in question. Then, the observable in equation 
(0) becomes a nonlinear function of the state, and may 
be termed a nonlinear nonclassicality witness (compare 



In the current Letter we restrict the class of investi- 
gated states g to those, for which P can be evaluated on 
an arbitrary real polynomial of x,y. The vector space of 
such polynomials will be denoted by R[a;,y]. Since any 
polynomial can be represented as a Fourier transform of 
appropriate sum of derivatives of the Dirac's delta func- 
tion, the sufficient condition for that is that P is a smooth 
function. We denote the space of such P's by V. 

Note, that since the test function / appearing in (@J) 
is entire, it can be almost uniformly approximated by a 
sequence of complex polynomials on C^. Hence |/r|^ = 
limjv^oo('«^ + vff) for some un,vn G K[a;,?/]. The al- 
most uniform convergence on the real plane allows us to 
interchange integration and taking the limit in Q [l^ . 
This leads to the main theorem of the present Letter: 

A state g with P E V is classical iff for every polyno- 
mial V 6 M.[x, y] 



d^aPv'^ = ti{Q : u(a,at)2 :} > 



(6) 



In fact the criterion © has already been used for a long 
time for detecting some important classes of quantum 
states. There are two examples of the application of ® 
known to the authors. The first one is the test for higher 
order quadrature squeezing [13 : g is squeezed to the 
order 2k if there exists a phase (j) £ [0, 2n) such that 



1 (2fc)! 

f^Jn\[2{k-i)V- 



\2{k-l) 



:) <0, 



(7) 



with := ae"'"* -|- ate'"^, AE^ := E^ ~ (E^) and the 
averages taken w.r.t. g. Obviously, {T)) has the form of a 
violation of © (we can always substitute v'^ there with 
finite sums of such terms) with the polynomial: 



W2k( 



1=0 



where 

d4x,y) ■.^2\x-{^) 



cosri 



y-(^)lsin0. (9) 



The witness W2k depends on the tested state g and hence 
is a nonlinear witness. 

The second example is the test for sub-Poissonian 
statistics of a^a (number squeezing): g is number 
squeezed if (: (Aa^a)^ :) < 0. The corresponding nonlin- 
ear witness is 



wp{x, y) := (x^ +y^ ~ (a^a))^ 



(10) 



Note that both nonlinear witnesses (jHl and (|l()|l are 
optimal in the sense that they are zero on the extreme 
points of A, as for any \a) all the moments of normally 
ordered deviations vanish. 
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From ©, we observe that for any v € R[a;,y], v"^ is 
positive semidefinite [psd), and so is every polynomial 
which is sum of such terms (we call such polynomials sos 
polynomials). One may ask if the converse is also true, 
i.e., if every psd polynomial is sos? 

This problem has been known in mathematics under 
the name of Hilbert's 17th problem. The answer is, quite 
surprisingly, negative: there are psd polynomials which 
are not sos ^] . For the case of 3 variables this happens 
for a degree m > 6. However, the explicit examples of 
psd, but not SOS polynomials are rare and were found 
quite lately. It is also worth mentioning that the con- 
nection between Hilbert's 17th problem and separability 
was established in 18]. 

In light of the theorem (jSJ, out of all psd polynomials, 
sos polynomials are sufficient to detect nonclassical states 
among the states with P € V. To illustrate how the 
theorem © works let us consider a specific example of 
sixth order Motzkin polynomial which is psd, but non- 
sos: 

M{x, y, z) := {x^ + - 3z^)x^y^ + . (11) 

Using a method origina ting from the witness techniques 
in entanglement study we construct a state g, de- 
tected by the polynomial M{x,y, ±1) [l^ . 

Out of the four zeros {(±1,±1)} of M{x,y,±l) we 
construct coherent states: ai :~ 1 -\- i, 012 '■— —1 + i, 
Q!3 :— Q!4 ■— a^. We pick the barycentric point, p, 
of the face J- := com/{S{a — ai), . . . , S{a — a^)} (conv 
stands for a convex hull) of the simplex A. Note, that 
the hyperplane defined by the witness : M(a, a\ ±1) : 

hM {PeP;J dxdy P{x, y)M{x, y, ±1) = 0} , (12) 

contains the face J- C A and hence the witness is optimal. 
Thus, to get the state detected by (|11|) . we mix g with a 
projector onto an arbitrary vector from it's range: 

. _ 4 

g:^—lj2\a,){a,\+e\i:){^\, (13) 

which for simplicity we choose to be: 

^:^l-{\a,) + \a*)), = 2[1 + c-\os{2)] . (14) 

Here < e < 1 and i G {1,2,3,4} is fixed, but the 
results presented below do not depend on it's particular 
value. Calculating the average of the polynomial l|ll|) 
using the expression we obtain: (: M{a,a^ ,±1) :) = 
(2/A^2)g-2pQg(2)e. Since cos(2) < 0, the state ^ is 
detected by M for e > 0. 

As a side remark, we note that the state (I13II is also de- 
tected by another example of psd, but non-sos polynomial 
- Choi-Lam polynomial 5'(a;, y, z) := x'^y'^ + y'^ z"^ + z'^x'^ — 
Sx^y^z^, as (: S{a, a^, ±1) :) = -{4/N^)e-^sm{2) e < . 



Before we explicitly construct a sos polynomial de- 
tecting (|13|l . let us first examine the physically relevant 
witnesses and (110(1 . A simple calculation gives that 
(: (Aa^a)^ :} > for any < e < 1. Examination of 
the witnesses 'W2k is more difficult and we have carried 
it out only numerically. We checked that up to the 14th 
order (k — 7) all the inequalities ^ are violated for any 
(j) & [0, 27r) and < e < 1 and hence fl^ is not squeezed 
up to the order of 14. Apart from that we used a modified 
version of W4: (x, y)^^^^ (a;, y)^ -I- Gd^g (x, ?/)^, depend- 
ing on three angles 4>i, 4>2 and (^3, also with no success. 
The question if g has even higher order squeezing is open. 

To construct a sos polynomial detecting g, note that 
M{x,y,±l) has only four zeros, and hence we can find 
a second order polynomial with the same zeros, which 
squared will give us the desired witness. Equivalently, we 
look for such a sos witness W that its hyperplane hw, 
defined as in ifT^ . contains T. We choose W{x,y) — 
(At^ + + Cxy + Dx + Ey + F)'^. The condition 
VF(±1,±1) ~ leads to a system of four linear equa- 
tions for A, . . . ,F. Its solution gives a family of witnesses 
WA,Bix,y) = {Ax^ +By^- A- B)^, where A^ + B'^^Q. 
The average of Wa,b in the state ifT^ is negative iff 
cos(2)((A + Bf - 4^2) -I- 4sin(2)A(A + B) < 0. As 
this equation possesses non-zero solutions, for example 
A = 0, B 7^ 0, the state g can be detected by a fourth 
order sos polynomial. 

This seems to be a generic feature, at least for the 
psd polynomials of degree m = 6. In this case from 
we know that if a psd polynomial has exactly ten zeros 
in PR^, than it cannot be sos. Fixing the variable z 
generally reduced the amount of zeros and hence permits 
to find a lower order sos polynomial with the same zeros. 

The methods described above, together with the cri- 
terion (jHJ, can be used to classify the states according 
to the degree of sos polynomial detecting them. Let us 
define a family of convex subsets of P: 

Sm-= f] {PeP;J d^aPw>0}, (15) 

where Em is the set of (inhomogeneous) sos polynomials 
of degree m. Theorem implies that A = p|j, S2k- It is 
also clear that S2 C E4 C • • • and hence S2 D S4 D ■ ■ ■ . 
We prove a stronger result: 

For any even m there exist nonclassical states detected 
by some witness from E„i, and not by any witness from 

^(m-2); that is ^2 2 ■5*4 2 ' ■ ' ■ 

Proof. Let us choose a generic w G E^. It has 
(r7i-|-l)(m-|-2)/2terms, as it is a sum of polynomials of de- 
gree < m. From the variety V{w) :— {{x, y);w{x, y) — 0} 
we pick n points (xi, yi), . . . , (x„, y„), m{m + l)/2 < 
n < {m + l){m + 2)/2, such that they do not lie on 
any variety of the lower order V{u), u G E(m_2)- We 
can find such points, as otherwise there would exist 
u e E(m-2), such that (xi, yi), . . . , (x„, y„) £ V{u). 



4 



However, with chosen n the latter condition leads to 
an overcomplete system of linear homogeneous equa- 
tions for the coefficients of u, which generically pos- 
sesses no solution. On the other hand, same condition 
for V{'w) yields under-determined system possessing non- 
trivial solution. Having such points we construct coher- 
ent states \xi + iyi), . . . , |x„ -I- iyn) and a face Tn G A 
spanned by them. For any g G 3-n we have then that 
tr{g : w(a, a^) :} = 0, whereas tr{^ : u(a, a^) :} > 
for all u € E(m_2)- Hence we can find such a con- 
vex combination g oi g and a projector onto some lin- 
ear combination of \xi + it/i), . . . , |a;„ + ij/n), such that 
trj^i : w{a^a^) :} < 0, while for all u G E(,„_2), 
iv{g : u{a,a^) :} > (from the continuity) . □ 

Summarizing, we have derived a family of classicality 
criteria of states of a quantum system, that require that 
P-representation averages of positive functions are pos- 
itive. For polynomial functions, we have related these 
criteria to 17-th Hilbert's problem: we have proven the 
theorem that all " generic" non-classical states (for which 
the P-representation averages of polynomials exist), can 
be detected by sos polynomials of sufficiently high de- 
gree; in this sense non-sos polynomials (whose existence 
was proven by Hilbert) are not necessary for classical- 
ity detection. We have also introduced the hierarchy of 
states implied by this theorem, and have introduced con- 
vex sets 52 2 'S'4 2 ■ • • of states detected by squares of 
polynomials of the 1st, 2nd,... order, and corresponding 
in this sense to decreasing degree of quantumness. 

We stress that our results have important experimental 
consequences. Our polynomial nonclassicality witnesses 
can be easily measured, allowing thus for direct detec- 
tion of quantumness and its degree for a given state. 
In this sense they are similar to entanglement witnesses 
that are nowadays commonly used for detection of en- 
tanglement . If one wants to check if a given state g 
is quantum, it is enough to measure normally ordered 
averages of squares of real polynomials of position q, 
and momentum p, or quadrature operators. In order to 
check the degree of quantumness (i.e. to check whether 
Q G S2k), one should determine normally ordered aver- 
ages of squares of real polynomials of the order k. Note, 
that for a given k this requires measurements of finite 
number of averages only. For instance, for k = 1 (squeez- 
ing), one needs to measure (q), (p), {: q^ :), (: :), and 
(: qp + pq :), and check if there exist A, B,C such that 
(: {Aq + Bp + C)^ :) < 0. For general k, one needs respec- 
tively fc(2fc -f- 3) measurements. Our results for the first 
time fully categorize states with respect to their degree 
of quantumness, and generalize the concepts of (higher 
order) squeezing or number squeezing as a signature of 
quantumness. 
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